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We use numerical and analytical results to construct a simple ansatz for the energy dynamical 
correlation function of the ID antiferromagnetic Heisenberg model. This is applied to compute the 
phonon assisted absorption spectra of magnetic excitations (spinons) in quasi-one dimensional spin 
1/2 insulators and show to reproduce very well recent infrared measurements in Sr2CuC>3. 
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After the solution by Bethe and many years of intense 
research, the spin 1 /2 Heisenberg chain is one of the most 
studied many-body problems. In spite of this little is 
known on its dynamical properties. This is because the 
extreme complexity of Bethe's wave functions makes it 
too hard to compute dynamical correlations. Much of 
the progress came from numerical studies which inspired 
an ansatz for the dynamical structure factor 0J^], a re- 
cent analytical evaluation of the exact two-spinon contri- 
bution H and low energy bosonization studies [||-|| . 

In this letter we use numerical and known analytical 
results to obtain a simple ansatz for the energy dynam- 
ical correlation function defined below. This correlation 
function is important when considering the problem of 
coupling with phonons || and in optical properties . 
As an application we compute the line shape for optical 
absorption of magnetic excitations assisted by phonons 
j^D and show how spinon physics appears naturally in 
the spectra. This is shown to reproduce very well recent 
measurements of the line shape by Tokura et al. PJTot in 
quasi-one dimensional cuprates (See Fig. ||). 

We will consider the 1-dimensional (ID) Hamiltonian, 
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is the bond-energy operator, A measures the anisotropy 
and we take periodic boundary conditions. We will con- 
sider A = (the XY model) and A = 1 (the antiferro- 
magnetic Heisenberg model). 

The energy correlation function is defined as 
G(p,ui) = ((SB p ; 5B_p))(u}) (in Zubarev's notation (TTI 
where 5B p is the Fourier transform of SBi = Bi — (Bi 



To warm-up we will compute G(uj,p) for A = 0. This 
is trivially done by using a Jordan- Wigner transforma- 
tion Jl2| to map the problem to a noninteracting ID free 
fermion gas. Excitations are free fermions and G(p, lo) 
becomes a Lindhard like response describing density like 
fluctuations. We get for the imaginary part, 

-2ImG XY (p,Lo > 0) = D XY (p, uj)M XY {p, w) (3) 

where following the approach of Ref. Q we defined a 
density of states D(p, ui) and a matrix element M(j>, ui) . 

Dxy [p, w) = — . = = (4) 
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with u>i (p) = Jsin(p), io* (p) = 2Jsin(p/2) and 
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Two-fold degenerate states are counted once in Eq. (Q) 
and this is taken into account by the factor of 2 in the 
matrix element (See Ref. [||). 

Notice that ImGxy(p,w) is only different from zero 



between to* Y and i 



,XY 



For a finite system the spectrum 



is given by a two-parameter set of discrete states with 
energy and pole strength given for p > by, 



Uxrip, q) = u£ (P) cos (<? ^ P/2) 
M XY (p,q) = -£jSm 2 (q~p/2) 



(6) 



where < q < p/2. Eq. ([|) corresponds to particle-hole 
excitations around the — ^ Fermi point. 

When we switch to the isotropic limit (A = 1) a simple 
analytic way of computing G(p, ui) is not available. Next 
we will follow Ref. O and use numerical results in finite 
chains, known low energy bosonization results and Bethe 
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ansatz results to construct a phenomenological ansatz for 
Im G(p, u)). 

In the isotropic limit there is no upper bound for the 
support of Im G(p, to) however the dominant contribution 
to lmG(p,uj) comes from singlet states between uj\{p) 
(the des Cloizeaux-Pearson [jl3| dispersion relation) and 
uj^ip) which are obtained from (p) and w£ Y (p) by 
changing J according to 

Jxy -> Jtt/2 (7) 

and we attached an XY subindex to avoid confusion lat- 
ter. This is analogous to the scenario found for the dy- 
namical structure factor [^) and corroborated by a recent 
analytical evaluation which retains only two-spinon ex- 
cited states ||. We illustrate the analogous behavior for 
the energy correlation function for a finite system in Fig. 

0- 




FIG. 2. The dots are the matrix elements (times N 2 ) vs. p 
from exact diagonalization of a 26 site Heisenberg ring. We 
had plotted the fit of N M(p, q) (the ansatz) as continuous 
functions of p for each allowed q in the 26-site chain. The 
intersection of the continuous curves gives the actual values 
for the finite chain. 




FIG. 1. Exact di2g%5aliMi8n 2 r'es 3 ults for ImG(p,w) in the 
uj — p plane for a 26 site lleisenberg ring. Excitations are 
given by the position of the circles and matrix elements by the 
radius. The dots are the XY energies with the renormalization 
of Eq. (Q). The lines are u>i(p) (lower) and u)2{p) (upper). 

In the thermodynamic limit these states form a two- 
parameter continuum degenerate with the triplet states 
that dominate the dynamical correlation function and 
have a density of states, D(p,u>) = ~* j7r ^ 2 (p, u>) [Q. 

To construct an ansatz for the energy correlation func- 
tion of the Heisenberg model we put 

-2lmG(p,u>Q)=D(p,u)M(p,u>). (8) 

Our ansatz for M(p, u) is given by, 

M(p,u) = M J x ^ J * /2 M(A^jL + B J — ) 

(9) 

For a finite system we use Eq 's (^) to construct M{p, q) 
in the same manner. In Fig. ^ we plot the matrix ele- 
ments of Fig. ^] as a function of p together with the best 
fit of M(jp, q) from the ansatz. The agreement is very 
good for all branches with significant spectral weight. 



It is interesting that one can fit all branches except 
the upper one with A ~ 3 and B — 0. The upper branch 
corresponds to the states closer to u>x (p) where the square 
root becomes important, i.e. the dominant correction to 
the XY matrix elements in most of the (p, u>) plane is 
given by the p independent factor. 

To determine A and B for an infinite system we have 
made a finite size scaling study and used the sum rule 
J2 P I duj(-2)ImG(p,uj) = In 2- (ln2) 2 to get A = 2.4, 
B = 0.6. 

Bosonization studies predict a divergence for 

— ImG(p — 7T, w) like w _1 which becomes (uj — cji(p))~ 1//2 
close but not at p =*k. The term proportional to B in 
M(j>,uj) ensures this behavior and is analogous to the 
form introduced by Muller et al. p| for the dynamic struc- 
ture factor . As mentioned before its importance here is 
quantitatively minor in most of the (p, uj) plane whereas 
in the case of the dynamic structure factor it was the 
dominant contribution. 

We next apply our results to the computation of the 
line shape for phonon assisted optical absorption of mag- 
netic excitations for a ID spin 1/2 isotropic Heisenberg 
chain. The absorption coefficient is given by |t],|8| , 

a(uj) — otobjl{bj — W||) , (10) 

where ao determines the oscillator strength and is given 
by 

ao = —=rj^ , (11) 

M is the O mass, Vc u is the volume associated with a 
Cu ion, uj\\ is the frequency of the Cu-0 stretching mode 
phonons, e± is the real part of the dielectric constant, 
<7a is an effective charge associated with this process and 
is of order \q&\ ~ e(3aJ/(SA 2 ), where (3 is the rate of 
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change of the Cu on-site energy due to a nearby O-ion 
displacement, a is the lattice constant, S — 1/2 is the 
magnitude of the spin and A the on-site energy difference 
between Cu and O. 

The function I is given by 

I{uj) = -16^sin 4 (p/2)ImG(p,w). (12) 
p 

I(u>) has a logarithmic singularity at to — Jtt/2 which will 
be convoyed to the optical absorption. The position of 
the singularity is a very robust feature which comes from 
the saddle point at p = 7r/2 in the exact des Cloizeaux- 
Pearson dispersion relation (See Fig. Q). The dis- 
tance between ui\\ and the singularity in the data pro- 
vides a quick way to obtain the microscopic parameter J 
directly from experiments. 

The two-dimensional version of this theory within 
an interacting spin-wave approach was recently applied 
to explain related mid-IR bands measured by Perkins et 
at. m spin 1/2 || and spin 1 jl5| insulating materials. 
Unlike estimations done in higher dimension J7|p||l^] the 
J obtained from the position of the singularity in ID 
is "exact" since the location of the singularity is known 
exactly. Here the only approximation is the Heisenberg 
model itself. 

In Fig. U we compare the calculated optical absorption 
line shape to recent optical data by Tokura and collab- 
orators P|l~C|] in Sr 2 Cu03. The agreement is excellent. 
We mention that the agreement, specially close to the 
singularity is not so good if B is set to zero. This suggest 
(in accordance with the numerical evidence) that the ex- 
ponent predicted by bosonization at low energies persist 
up to energies quite far from its range of applicability. 

J = 0.246eV was adjusted and found to compare very 
well with the value obtained in susceptibility measure- 
ments 0. 

The value of «o which was fitted to the observed spec- 
tral weight, can be used to estimate the effective charge 
oa- We find |oa| — 0.4e whereas in the 2D cuprates 
\q^A ~ O.le was found. This is not surprising since the 
processes that give rise to qa are quite similar to the ones 
that give rise to J and the latter is found to be a factor 
of 2 larger in the ID compounds. 

Both the singularity and the cutoff at W2(tt)+ re- 
flect the dominance of the spectra by two-spinon states 
analogous to the two-spinon states considered in Ref. ||] 
for the dynamical structure factor. 

In conclusion we presented a simple ansatz for the en- 
ergy dynamical correlation function of the Heisenberg 
model. The ansatz reproduces very well numerical data 
and the low energy behavior predicted by bosonization 
studies and provides new insight on the dynamical prop- 
erties of the ID Heisenberg 
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FIG. 3. Experimental absorption coefficient from Ref. || 
(solid line) and theoretical line shape (dashed line) in 
Sr 2 Cu0 3 . Parameters are J = 0.246 eV, wy = 0.08 eV and 
ao = 132cm -1 . aj|| was taken similar to the value in the 
2D materials. A linear background was subtracted from the 
experimental data. 

After this work was completed we became aware of 
Ref. JkJ which contains the experimental data and a the- 
oretical analysis of the line shape qualitatively similar to 
ours. They predicted the same position for the singu- 
larity however their matrix elements where computed in 
the XY limit. Although at first glance the theoretical 
line shape resembles the experiment they have missed 
the factor u in Eq. (|l^) and the p-dependent form factor 
in Eq. (0). When they are included the resemblance is 
partially lost. Additionally they get a cusp singularity 
instead of a logarithmic singularity. 
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